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Spectral singularities that spoil the completeness of Bloch-Floquet states may occur in non- 
Hcrmitian Hamiltonians with complex periodic potentials. Here an equivalence is established be- 
tween spectral singularities in complex crystals and secularities that arise in Bragg diffraction pat- 
terns. Signatures of spectral singularities in a scattering process with wave packets are elucidated 
for a 'PT-symmetric complex crystal. 
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I. INTRODUCTION 



The Dirac-von Neumann formulation of quantum me- 
chanics prescribes that the Hamiltonian H of a physical 
system must be Hermitian. This requirement ensures a 
real-valued energy spectrum and a unitary (probability- 
preserving) temporal evolution. In the last decade, 
complex extensions of quantum mechanics which relax 
the Hermiticity constraint have been proposed (see e.g. 
[l|, [|j]). Indeed, a diagonalizable Non-Hermitian Hamil- 
tonian (NHH) having a real spectrum is enough to con- 
struct a unitary quantum system, provided that the inner 
product of the Hilbert space is properly modified @ • In 
particular, Bender and collaborators showed that a NHH 
possessing parity-time (VT) symmetry may serve to de- 
velop a complex extension of quantum mechanics below 
a phase transition (symmetry- breaking) point [H, 0| • Ex- 
perimental realizations of VT Hamiltonians have been 
recently proposed in optical media with a complex re- 
fractive index 0411; and the first observation of VT- 
symmetry breaking has been reported in a passive optical 
waveguide coupler jH. 

For a NHH, the reality of the spectrum does not ensure 
diagonalizability, which may be prevented by the pres- 
ence of exceptional points in the point-spectrum of H 
[loj . or of spectral singularities in the continuous part 
of the spectrum [Ill4l3| |. Exceptional points, also re- 
ferred to as Hermitian degeneracies (see |10|j) , correspond 
to degeneracies where both eigenvalues and eigenvectors 
coalesce as a system parameter is varied. Other singu- 
lar points may also occur in low-dimensional (e.g. ma- 
trix) non-Hcrmitian Hamiltonians, such as branch points 
in the complex plane in which eigenvalue degeneracy do 
not correspond to lack of completeness of the spectrum 
[l4j . Such singularities have various physical implica- 
tions, which have been investigated in different physi- 
cal fields (see for instance @, 0, OH [H| and references 
therein). They play an important role in the study of 
open quantum systems, particularly in relation with the 
resonance states (for a recent review see [16]). But they 
must not be confused with spectral singularities as de- 
fined in 

una. 

Unlike exceptional and branch points 
that can be present for non-Hermitian operators with 
a discrete spectrum, spectral singularities are exclusive 



features of certain non-Hcrmitian operators having a 
continuous part in their spectrum [13| . Following jlll - 
[l3j . spectral singularities refer to divergences (poles) of 
the resolvent on the continuous spectrum of H, which 
do not correspond to square integrable eigenfunctions. 
The physical implications of such spectral singularities 
have been investigated so far for wave scattering from 
complex potential barriers [H, EH E3, EH and shown 
-notably by Mostafazadeh- to correspond to resonance 
states with vanishing spectral width [l?], [Hi- Another 
physically relevant class of NHHs with continuous spec- 
trum is provided by complex periodic potentials (22l - |24j . 
Complex crystals have been investigated in different ar- 
eas ofphysics, ranging from matter waves jl9H2l| to op- 
tics [a S US [HI, HH ■ As compared to ordinary crys- 
tals, complex crystals exhibit some unique properties, 
such as violation of the Friedel's law of Bragg scattering, 
double refraction, nonreciprocal diffraction, and anoma- 
lous transport [(| [l^ - [2ll [25| , which make them a rather 
unique class of synthetic materials. Spectral singularities 
for periodic non-self-adjoint Schrodinger operators have 
been studied by mathematicians [27| , however their phys- 
ical implications in connection to complex crystals have 
not been explored yet. 

It is the aim of this work to investigate the onset of spec- 
tral singularities in complex crystals and to show their 
physical impact in Bragg scattering processes. In par- 
ticular, we prove rather generally that spectral singular- 
ities are associated to a secular growth of plane waves 
diffracted off the crystal when the incident angle is an 
integer multiple of the Bragg angle. With reference to a 
specific PT-symmetric complex crystal, previously con- 
sidered in Ref.[6j, we show that spectral singularities 
occur at the VT symmetry breaking point and can be 
revealed in a diffraction experiment using a spatially- 
confined wave packet that excites the crystal at normal 
incidence. Unlike plane wave excitation, broadening of 
the angular spectrum for a wave packet is shown to lead 
to a saturation of the secular growth of the wave ampli- 
tude. 
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II. SPECTRAL SINGULARITIES IN COMPLEX 
CRYSTALS AND BRAGG SCATTERING 

A. Bloch-Floquet states and spectral singularities 

Wave dynamics in a complex one-dimensional crystal 
is governed by a Schrodinger-type equation, which in di- 
mcnsionless form can be written as 

idti>{x, t) = -d%ip + V(x)il> = HiP (1) 

where V{x + a) = V(x) is the complex potential of pe- 
riod a. Physically, Eq.(l) describes diffraction of matter 
or optical waves by a complex lattice [see Fig. 1(a)], where 
t is a fictitious time related to the propagation distance 
inside the crystal [f| [IH, l25l] . The geometry of the 
spectrum of H is studied e.g. in [24j ; here we will assume 
an entirely real energy spectrum, however at this stage no 
requirement about VT symmetry of H is needed. Like in 
ordinary crystals, from the canonical form of the transla- 
tion operator [26( it follows that any eigenfunction <fi(x) 
of H (H<f) = E<f>) which is bounded at x — > ±00 is a 
linear combination of Bloch-Floquet type solutions sat- 
isfying the condition <f)(x + a) — <f>(x) exp(iqa), where 
— 7r/a<g<7r/aisan arbitrary real number that varies 
in the first Brillouin zone. The set of Bloch-Floquet 
eigenfunctions with energy E a (q) can be expanded in se- 
ries of plane wave basis \q, n) = (27r) -1 / 2 exp[i(q+nkB)x] 
as 

Mx, q ) = J2^nHq)\q,n), (2) 

n 

where ks = 27r/a is the Bragg wave number, a is the 
band index, and n = 0,±1,±2,.... In Eq.(2), w (a \q) = 

{w^"\q)} are the eigenvectors, with corresponding eigen- 
values E a (q), of the matrix %{q) defined by 

Un.miq) = (q + nk B ) 2 S ni m + V n -m (3) 

with V{x) = J2n exp(jfcs^). Let J{q) be the Jor- 
dan canonical form of 'H{q) 1 to which Tl{q) can be 
reduced by a similarity transformation, i.e. T-L{q) = 
T{q)J{q)T^ 1 {q). In the band computation, degenerate 
eigenvalues of H(q) are counted by their geometric (and 
not algebraic) multiplicity, i.e. by the number of Jor- 
dan blocks that represent the eigenvalue, and the eigen- 
vectors w( a (q) are thus linearly independent. Note that 
the matrix representation of H in the plane-wave basis 
is given by (q',n'\Hq,n) = H n > , n {q)5(q - <?')• Like in 
ordinary crystals, the eigenvalues can be ordered such 
that E a {—q) = E a (q), and for any fixed value of <?, with 
q 7^ 0, — 7r/a, the energies E a {q) are all distinct. At q = 
or at q — —it /a, eigenvalue degeneracy, with algebraic 
multiplicity not larger than 2, is allowed (see e.g. [26|). 
For the problem of spectral singularities, a key role is 
played by defective eigenvalues of H(q), at either q = 
or q = — ■n/a, i.e. eigenvalues whose geometric multi- 
plicity is smaller than their algebraic multiplicity. More 



precisely, we will prove that the completeness of a peri- 
odic NHH admitting a purely continuous spectrum, i.e. 
the absence of spectral singularities in complex crystals, 
is equivalent to the issue of completeness (i.e. absence 
of defective eigenvalues) for the non-Hermitian matrix 
Hamiltonian %{q) at q = and q = —n/a. Note that 
the discrete NHH problem defined by the matrix H(q) 
at q = 0,— it /a can be obtained from the original prob- 
lem (1) provided that the functional space is restricted to 
the space of functions satisfying the periodic boundary 
conditions ip(x + a,t) = ±ip(x, t). Moreover, the appear- 
ance of a defective eigenvalue of H(q) (when a control 
parameter is varied) corresponds to the occurrence of an 
exceptional point in the discrete spectrum of H(q). The 
correspondence between spectral singularities in the con- 
tinuous spectrum of H and exceptional points (defective 
eigenvalues) in the point spectrum of the matrix H(q) is 
more precisely established by following theorem: 
Theorem 1. The set of Bloch-Floquet eigenfunctions (2) 
is complete, i.e. if is diagonalizable, if and only if the 
matrix 'H{q) does not have defective eigenvalues at q = 
or q = — 7r/a. In other words, spectral singularities in 
the continuous spectrum of H correspond to exceptional 
points in the point spectrum of the matrix T-L{q) at the 
center (q — 0) or at the edge (q = — it /a) of the Brillouin 
zone. 

Proof. Let us first prove that {0 Q (x, q)} is a complete 
set of (improper) functions if and only if T-L{q) is diag- 
onalizable for any arbitrary value of q inside the first 
Brillouin zone. Suppose that {4> a (x,q)} is a complete 
set. Then for any q in the first Brillouin zone and 
for any integer n, there exists a set of complex num- 
bers {c n , a (q)} such that \q,n) = J2 a c n , a (q)(f> a (x, q), i.e. 

\q^ n ) = Em (j2 a c n,a(q)wm\qfj \q,m), where we used 

Eq.(2). Hence Y Ja c n,a{q)wm\q) = S n ,m- Since the pre- 
vious relation holds for any integer value n, the set of 
eigenvectors {w( Q )(g)} is complete, i.e. T-l(q) is diagonal- 
izable. Suppose now that 'H(q) is diagonalizable, and let 
us show that {cj) a (q, x)} is a complete set. In fact, for 
an arbitrarily assigned square integrable function f(x), 
let F(k) be the Fourier transform of f(x) and F n (q) — 
F(q + nkB), where q varies in the first Brillouin zone and 
n = 0,±1,±2,.... Then f{x) = Y. n S%%dqF n {q)\q,n). 

Since V.(q) is diagonalizable, the eigenvectors w( Q )(q) 
form a complete set and one can thus determine a set 
of functions c a (q) such that F n (q) = ^ a c a (q)v}k\q). 
Hence f(x) = J2 a !-^ a d( 7 c «(<?) E„ w< n ] (q)\q,n) = 
E Q Jli a / a dqc a (q)(p a {x,q), i.e. f(x) can be decomposed 
as a superposition of Bloch-Floquet eigenfunctions. Since 
f(x) is arbitrary, it follows that {4> a {x, q)} is a complete 
set. The theorem is finally proved after observing that, 
for q =/= 0, — 7r/a, the eigenvalues E a (q) are distinct, and 
thus H(q) is diagonalizable. 

Similarly to spectral singularities found in scatter- 
ing complex potentials (T3 |. a physically relevant conse- 
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quence of spectral singularities is to prevent the construc- 
tion of a biorthogonal eigensystem for H, i.e. to resolve 
the identity in terms of the biorthogonal basis associated 
to H. By definition, H is diagonalizable if cj> a (x,q), to- 
gether with a set of (generalized) eigenfunctions <j>a(x, q) 
of the adjoint H\ form a complete biorthogonal sys- 
tem, i.e. they satisfy {4> a (x, q)\<f>^(x, q')) = S a ^5(q - q') 

and E a j:i a /a d q \0 a (x >q ))(<f>i(x 7q )\ = X. As fft = 
— d"l + V*(x), it can be easily shown that one has 
4>i( x , 0) = N a {q)<t>* a {x, ~q)> where J\f a (q) is a multiplying 
term that needs to be determined. A direct computation 
of the scalar product ((f> a (x, q)\(jyg(x, q')) using expansion 

(2) yields (<f> Q (x,q)\<f>l(x,q')) = 6{q- q')S a , p Af a (q)V a (q), 
where we have set T> a (q) — (w^ (q)\w^ a ^ (q)) and 
w ( Q )t(q) = W ( Q )*(— qr). Therefore, completeness of 
the biorthogonal system is ensured by letting Af a (q) — 
1/V a (q) provided that that V a (q) ^ 0. Because V a {q) 
vanishes if and only if the energy E a {q) is a defective 
eigenvalue of the matrix %{q) 1 from theorem 1 one con- 
cludes that a spectral singularity prevents the construc- 
tion of a biorthogonal eigensystem for H. 
Finally, it is also worth mentioning the following the- 
orem, which shows that the defective eigenvalues of H 
correspond to the spectral singularities of H, defined as 
the divergence points of the resolvent G(z) = (z — 
on the continuous spectrum. 

Theorem 2. Any defective eigenvalue E of % is a diver- 
gence point for the resolvent G{z) — (z — If) -1 . 
Proof. Let |x) and \<p) be two square integrable functions 
of the Hilbert space, and let us introduce the complex 
function of z, G x>{p (z), defined by G XiV (z) = (x\G(z)ip). 
To prove the theorem, it is enough to show that there 
exists at least a couple of functions |x) and \<p) such that 
G x>t p{z) is unbounded as z —¥ Eq, Eq being a defective 
eigenvalue of %. To this aim, after expanding |x) and 
\tp) on the plane- wave basis \q,n), one can readily show 
that 

/7T j a 
%ffl(?)^(?)Km, n (2,9)> (4) 

where Xm(o) = (q, m \x), fn(q) = (q, n \f) and 
K(z,q) = [Z - H(q)}- 1 = T(q)[z - J(q)]- l T- l {q) 
is the resolvent of the matrix H(q). Note that the 
spectral functions Xm(q) and <p n (q) are simply re- 
lated to the Fourier transforms of |x) and \tp) as 
described in the proof of Theorem 1, and hence the 
inversion relations 1x0*0) = J2n 1-1% d 1Xn(q)\q, n) 
and \<p(x)) = J2nI-i'/ a d 1 ( Pn(q)\q,n) hold. If E Q is 
a defective eigenvalue of H(q), say at q = 0, with 
algebraic multiplicity 2 and geometric multiplicity 1, 
then there exists one element of the matrix lZ(z,q), say 
T^"m ,n ( z ) q)i which has a second-order pole at z = Eq 
when q — 0. This follows from the well-known form of 
the resolvent [z — J^(q)]^ 1 of a Jordan matrix possessing 
a Jordan block with dimension 2. In the neighbor- 



hoods of q = and z — E , 72. m0i „ (z, q) behaves like 
~ l/[(z — Eq — aq)(z — Eq + aq)], where a is a constant 
and Eq ± aq are the two eigenvalues that cross and be- 
come a defective eigenvalue at q = 0. Let us now choose 
the spectral functions Xm{q) to vanish for m mo, 
and <fn(q) to vanish for n ^ tiq; conversely, Xm (q) and 
ip no (q) are assumed to be nonvanishing and bounded 
functions of q in the interval (— 7r/a, tt/o). The functions 
|x) and \ip) in direct space that yield such spectral 
functions are determined by the corresponding inversion 
relations given above. For instance, if Xm (?) and ip no (q) 
assume a constant value in the interval (— 7r/a, 7r/a), 
then 1x0*0) — A/"[sin(7rx/a)/(x)] exp(zfcsmox) and 
|<^(a;)) = J\f[sm(irx/a)/(x)]exp(ikBnox), where Af is a 
normalization constant. Under such a choice, the sum 
in Eq.(4) reduces to the only term m = tuq and n = uq. 
Owing to the behavior of Tt(z, q) near z = Eq and 
q = 0, the corresponding integral on the right-hand-side 
of Eq.(4) diverges as z — > Eq ± i0 + , i.e. G XiV (z) is 
unbounded in the neighborood of z = Eq, which proves 
the theorem. 



B. Spectral singularities and secular Bragg 
diffraction 

An important physical implication of spectral sin- 
gularities in complex crystals is the appearance of a 
secular growth of the amplitudes of waves scattered 
off the lattice when it is excited by a plane wave at 
special incident angles. Such an anomalous behavior 
has been previously noticed by Berry [2(| [2l| for certain 
absorptive potentials; here we prove that this is a 
very general feature of complex crystals related to the 
existence of spectral singularities. Namely: 
Theorem 3. Let ip(x,t) = cxp(—iHt)ipQ(x) be the Bragg 
diffraction pattern corresponding to crystal excitation, 
at t = 0, with a plane wave VK^jO) = exp(ifcx) of wave 
number k. Then H has spectral singularities if and only 
if for some integer n and k = nfc^/2 the solution ij)(x,t) 
contains secular (linearly growing) terms in t. 
Proof. Let us set k = q + lgkB, where Iq is an integer and 
— n/a < q < n/a. Then the solution ijj{x,t) of Eq.(l) 
with the initial condition ip(x,0) = cxp(ikx) is given by 
ip(x,t) = Ylil-oo c i(*) ex P[*(<? + lks)x], where cj(t) sat- 
isfy the coupled equations i(dci/dt) = ^ m T-Li,m(q)c m {t) 
with the initial conditions c;(0) = 8i t i . Physically, the 
coefficients q are the amplitudes of diffracted waves 
at various orders at plane t [see Fig. 1(a)]. Hence 
ci(t) = Mi,i (q,t), where M(q,t) = exp[—itH(q)]. The 
exponential matrix M. can be calculated from the Jordan 
decomposition of V. as M = Tcxp(— itJ')T~ 1 . Having 
in mind the form of the exponential of a Jordan matrix 
pjj , it follows that secular growing terms in some of the 
coefficients ci(t) [and hence in ip(x, t)] appear if and only 
if H(q) has at least one defective eigenvalue, the largest 
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growing term being ~ t p where p is the (maximal) dif- 
ference between the algebraic and geometric multiplicity 
of defective eigenvalues. In our case, p — 1 because the 
maximal algebraic multiplicity of any eigenvalue of H(q) 
is 2; moreover, since may have defective eigenvalues 
only for q = or q = —ir/a, secular growing terms may 
appear solely when the wave number k of the exciting 
plane wave is an integer multiple of fcs/2. The absence 
of secular growing terms in q (t) for any excitation wave 
number k = nk B j2 (n = 0,±1,±2,...) implies that the 
matrix J is diagonal, which ensures the lack of defective 
eigenvalues of T~L and thus of spectral singularities of H 
according to Theorem 1. 




FIG. 1: (color online) (a) Schematic of Bragg scattering of a 
plane wave off a complex crystal, (b) Band structure of the 
the complex crystal V(x) = Vo exp(jfc_ex). The circles mark 
the spectral singularities inside the continuous spectrum. 



III. BRAGG SCATTERING IN VT COMPLEX 
CRYSTALS AND WAVE PACKET DYNAMICS 

Let us specialize the previous results to the case of a 
■pT-symmetric lattice, for which V(—x) = V*(x). Let 
Vr(x) and AV/(x) be the real and imaginary parts of the 
potential, respectively, where A > measures the anti- 
Hermitian strength of H. The spectrum of H is real 
for A < A c , where A c > defines the symmetry break- 
ing point. According to the previous analysis, complex- 
conjugate pairs of eigenvalues for H(q) should appear as A 
is increased from below to above A c . The typical scenario 
that describes symmetry breaking in a finite-dimensional 
VT matrix is the appearance of an exceptional point via 
the merging of two real eigenvalues into a single real and 
defective eigenvalue at A = A c (see, for instance, [71. |29|). 
We note that the existence of such a branching point for 
a certain class of non-hermitian matrices as a control pa- 
rameter is varied wasproven in a rather generally way 
in Ref.(30| (see also 0]). We may thus conjecture that 
for a VT complex crystal symmetry-breaking is accompa- 
nied by the appearance of spectral singularity, which arise 
from defective eigenvalues of % at g = or g = — w/a. 
This scenario is in agreement with numerical or analyt- 
ical results obtained from band computation of specific 
complex periodic potentials (see, for instance, 0, Q). 

As an example, let us consider the VT crystal defined 
by 

Vr{x) = V cos(27ra;/a) , Vi(x) = V sin(2vra;/a), (5) 

which has been recently considered to highlight unusual 
diffraction and transport properties of complex optical 
lattices @, [25[. In this case, A c = 1 @ and at the sym- 
metry breaking point one has V(x) = Vq exp(iknx\ a pq- 
tential which is amenable for an analytical study [20l . |22| . 
For this potential, H has a block diagonal form, namely 
T-L n , m = (n + k B q) 2 S n ^ n + Vo5 m , n +i, and its eigenval- 
ues are simply the elements on the main diagonal, i.e. 
E a (q) = (q + ak B ) 2 (a = 1 ±L±2,...). This means 
that, as previously noticed [6, 22, 25], the band structure 
of the VT potential V = Voexp(ikB%) coincides with 



the free-particle energy dispersion curve E = k 2 , period- 
ically folded inside the first Brillouin zone [see Fig. 1(b)]. 
The eigenvalues are distinct for q ^ 0, — ir/a. At the 
crossings of the folded parabolas of Fig. 1(b), i.e. at 
q = and at q — —n/a, one has E_ a {q) = E a (q) and 
E\- a {q) — E a (q), respectively, i.e. the eigenvalues coa- 
lesce in pairs and become defective. Therefore, the con- 
tinuous spectrum of H, E > 0, contains a sequence of 
spectral singularities at E n — (n/cs/2) 2 , n = 1,2,3,... 
[see Fig. 1(b)], which spoil the completeness of the Bloch- 
Floquct cigcnfunctions. The defective nature of degen- 
erate eigenvalues at q = and q = —n/a can be readily 
proven by direct calculation of the eigenvectors w*-"' of 
%. Note that, as wave scattering from complex potential 
barriers enables a finite number of spectral singularities 
in the continuous spectrum (l2l . Il7| , in our example the 
number of spectral singularities is countable but infinite. 
According to Theorem 3, a secular growth of Bragg 
diffraction pattern for a plane wave that excites the crys- 
tal at normal incidence (or tilted by an angle which is an 
integer multiple of the Bragg angle) provides a distinctive 
signature of the appearance of spectral singularities at 
the VT symmetry-breaking transition point A = A c = 1 . 
However, in any experimental setting aimed to observe 
such a secular growth, the wave that excites the crystal 
is always spatially limited or truncated, and it is thus of 
major relevance to investigate the impact of spectral sin- 
gularities on the evolution of a wave packet with a broad- 
ened angular spectrum, an issue which was not consid- 
ered in previous works by Berry. Here we investigate the 
Bragg diffraction of a wave packet with a broadened mo- 
mentum distribution, ijj{x, 0) = f dkF(k) exp(ikx), and 
show that spectral broadening leads to a saturation of 
the secular growth of scattered waves. Such a saturation 
behavior is basically due to the fact that spectral sin- 
gularities are of measure zero (they are a countable set 
of points embedded in the continuous energy spectrum 
E > 0). For the sake of simplicity, we consider a shallow 
lattice and a wave packet with a narrow spectrum F(k) 
centered at k = — fcs/2 of width Ak <§; k B ■ Following the 
same lines detailed in the proof of Theorem 3, one can 
show that the diffraction pattern ip(x,t) can be written 
as the interference of different wave packets describing 
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one-side diffraction at various orders, namely 
ip(x, t) = ^ {x, t) + VqMx, t) + v$M*> *) 



(6) 



where ip n (x,t) — J dkF(k)c n {k,t) exp[i(fc + nks)x\ 
with Co(fc,i) = exp(— ik 2 t) and c n (k,t) = 
—i Jq d£cn-i(k,£) exp[i(k + nkB) 2 {^ — t)} for n > 1. For a 
shallow lattice, i.e. for \Vq\ <C 1, we can limit to consider 
the first two terms on the right hand side of Eq.(6). The 
leading term, ipo(x,t) — J dkF(k) exp(—ik 2 t + ikx), is 
simply the freely-diffracting wave packet that one would 
observe in the absence of the crystal and that propagates 
with a constant speed v = dx/dt = kg. The expression 
of ipi(x,t) is more involved, however its asymptotic 
behavior for t 3> jr/(fcsAfc) assumes a rather simple and 
physically interesting form, namely 



m I ks\ ( knx k n t \ I x 
^ ~ --F ( -- £ ) exp ( i-^- - i^L ) $ ( — 



ks 

where $(£) 




1 for |£| < 1 and <I>(0 
(b) 



(7) 
for 
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FIG. 2: (color online) Evolution of a Gaussian wave packet 
in the complex lattice V(x) = Vo exp(ikBx) for Vb = 0.2, 
a = 1 and w = 80. In (a) a snapshot of |^(a;,t)| is shown, 
whereas in (b) the profiles of [^/>(a;,t)| at a few values of t are 
reported. In (b), the dashed lines correspond to the wave 
packet evolution in absence of the lattice, whereas the dotted 
curves correspond to |^i(a;,i)| as predicted by Eq.(7). 

|£| > 1. Equation (7) shows that, owing to wave packet 
broadening in momentum space, the secular growth 
with t of the diffracted beam saturates to the value 
~ (n/kB)\F(— fcs/2)| while the beam assumes a square 
shape whose width w s spreads in space with a con- 
stant speed v = dw s /dt — ks- Note that v equals the 
translation speed of the freely-diffracting wave packet 
ipo ■ This behavior is confirmed by direct numerical 
simulations of Eq.(f), as shown in Fig. 2. The figure 
depicts the evolution of |?/>(:r,i)| for a Gaussian wave 
packet ij}{x,0) — exp[— (x/w) 2 — ikBx/2] with spectrum 
F(k) = [w/(2 v / ^)]exp[-(fc + fc B /2) 2 w 2 /4] that excites 
the crystal at t = 0. The evolution of the freely diffract- 
ing wave packet | tpo \ and of the asymptotic behavior of 
\tpi\ predicted by Eq.(6) are also shown for comparison. 
Note the formation of interference fringes on the left side 
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FIG. 3: (color online) Saturation of secular Bragg scattering 
for a Gaussian wave packet that excites the complex lattice of 
Fig. 2 at normal incidence. The curves show the evolution of 
the maximum wave packet amplitude ip m (t) versus t for a few 
values of input wave packet spot size w. Curve 1: w = 40; 
curve 2: w = 80; curve 3: w — 150; curve 4: w = 300. 



of the wave packet, which arise from the interference of 
tjjo and ipi. The spectral-broadening- induced saturation 
of the secular growth of the wave packet is clearly shown 
in Fig. 3, where the behavior of the maximum amplitude 
ipm(t) of the wave packet versus t, defined by 

i>m(t) = maxv\ip(x,i)\, (8) 

is depicted for a few decreasing values of the wave 

(a) 

160 




FIG. 4: (color online) (a) Evolution of a Gaussian wave packet 
[modulus of ip{x,t)] in the complex lattice (5) for Vb = 0.2, 
a = 1, w = 80, and below the VT symmetry breaking point 
(A = 0.9). In (b), the corresponding evolution of the maxi- 
mum wave packet amplitude %p m (t) is shown. 

packet input spot size w. Note that in the early stage of 
the dynamics the peak amplitude linearly increases with 
t, as expected for a plane wave according to Theorem 
3. The linear growth then saturates to a steady state 
value and, as Fig. 3 clearly shows, the saturation process 
occurs earlier for wave packets with smaller input size 
w. It should be pointed out that saturation of ij) m (t) 
to a steady-state value, as predicted by Eq.(7) and 
confirmed by numerical simulations depicted in Figs. 2 
and 3, is indeed a signature of spectral singularities of 
the underlying Hamiltonian that arise at A = A c = 1. 
This is clearly shown in Figs. 4 and 5, where a typical 
wave packet evolution is reported for the complex 
crystal defined by Eq.(5) either below (A = 0.9A C = 0.9, 
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FIG. 5: (color online) Same as Fig. 4, but above the VT 
symmetry-breaking point (A = 1.1). 



point, which are revealed as the saturation of the 
wave packet amplitude growth to a steady-state (non- 
decaying) value. If the lattice is realized in a periodic 
dielectric medium as discussed in Q , at light wavelength 
A = 1.5 /im, assuming a bulk refractive index n s = 1.5 
and a maximum refractive index change (both real and 
imaginary parts) of An = 2 x 10~ 4 , the simulations of 
Fig. 2 correspond, as an example, to an optical lattice 
with spatial period a ~ 6.2/im excited by a Gaussian 
beam of size w ~ 493 fiia; the spatial units along the x 
and t axes are l x ~ 6.2 (im and It — 477 /im, respectively. 
The same scales hold for the simulations shown in Figs. 4 
and 5. 



Fig.4) and above (A = 1.1A C = 1.1, Fig.5) the VT 
symmetry breaking transition point. In both cases, H is 
diagonalizable and there are not spectral singularities. 
Note that below the VT symmetry breaking point, 
the energy spectrum is real-valued and, after an initial 
increase, ipm(t) does not settle down to a steady-state 
value, rather it tends to monotonically decay (Fig.4) like 
in an ordinary (Hermitian) crystal. On the other hand, 
above the VT symmetry-breaking point (A = 1.1), the 
energy spectrum contains pairs of complex-conjugate 
eigenvalues, and i> m (t) monotonically increases, as 
shown in Fig.5. Therefore, the scattering process with 
wave packets can be used to highlight the appearance 
of spectral singularities at the VT symmetric breaking 



IV. CONCLUSIONS 

In this work, it has been shown that Bragg scatter- 
ing in complex crystals provides a physically important 
process to visualize the appearance of spectral singular- 
ities in non-Hermitian Hamiltonians with complex peri- 
odic potentials. In particular, clear signatures of spectral 
singularities could be gained in a diffraction experiment 
with wave packets. These results may suggest the in- 
vestigation of possible experimental systems to observe 
spectral singularities in complex lattices, such as pho- 
tonic systems [y, l9( . 
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